Structure, bonding and magnetism in cobalt clusters 
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The structural, electronic and magnetic properties of Con clusters (n =2—20) have been investi- 
gated using density functional theory within the pseudopotential plane wave method. An unusual 
hexagonal growth pattern has been observed in the intermediate size range, n =15—20. The cobalt 
atoms are ferromagnetically ordered and the calculated magnetic moments are found to be higher 
than that of corresponding hep bulk value, which are in good agreement with the recent Stern- 
Gerlach experiments. The average coordination number is found to dominate over the average bond 
length to determine the effective hybridization and consequently the cluster magnetic moment. 

PACS numbers: 75.75.+a, 36.40.Cg, 61.46.Bc, 73.22.-f 



I. INTRODUCTION 

Study of finite size clusters is an important means of 
understanding how magnetic behavior evolves in reduced 
dimensionality. The 3d transition metal (TM) elements 
are characterized by their unfilled d-shell, which gives 
rise to their magnetism and many other interesting phys- 
ical and chemical properties. Early transition metals are 
nonmagnetic in bulk solids, and only Fe, Co and Ni are 
known to be ferromagnetic among the 3d metals. How- 
ever, the small clusters of all the early transition metals 
are magnetic and those of late transition metals possess 
magnetic moments enhanced from their bulk values due 
to their spatial confinement. Local spin density based 
calculation^ showed that the face centered cubic (fee) 
phase is the lowest energy state for bulk Co in the param- 
agnetic phase, whereas the magnetic order stabilizes the 
hexagonal close packed (hep) phase as the ground state. 
This indicates a strong correlation between the stable 
structure and magnetism. Although, a metastable anti- 
ferromagnetic state exists for bulk Co, the ferromagnetic 
state is found to be the most stable for all crystal struc- 
tures. This is unlike other 3d transition metals Cr, Mn, 
and Fe which have a stable antiferromagnetic structure 
in their fee phase. This means that not only the crystal 
structure, but also the electronic configuration controls 
magnetism. In the present communication, we focus on 
Co„ clusters to understand this interplay. 

The magnetic properties of bare Co„ clusters were 
first investigated via Stern-Gerlach (SG) molecular beam 
deflection experiment by Bloomfield and co-workers for 
C020— C0215 clustersSii^ and by de Heer and co-workers 
for C030— C0300 clusters j^i^i^ These studies showed that 
in the temperature range of 77—300 K, the Co„ clus- 
ters display high-field deflections, which are character- 
istic of superparamagnetic behavior. The superparam- 
agnetic model for free clusters was recently revisited by 
Xu et alJ- and they proposed that adiabatic magnetiza- 



tion together with avoided Zeeman levels crossing in iso- 
lated clusters can lead to the same high-field beam de- 
flection behavior as observed in the superparamagnetic 
spin relaxation. However, both the models predict the 
same high temperature limiting form for magnetization 
as given by the Curie law^i^ (M) = ^"^B/ikT, where ^ 
is the cluster magnetic moment, B is the magnetic field 
and T is the cluster temperature. The intrinsic per-atom 
magnetic moment for small Co„ clusters was found to 
be substantially larger than the bulk value^'^i^i^'^ii'^ and 
generally decreases with increasing cluster size, eventu- 
ally reaching the bulk value at ~ 500 atomsi^ The en- 
hancement in the magnetic moment in small clusters has 
been attributed to the lower coordination of the surface 
atoms resulting in a narrowing of the d-bands and hence 
greater spin polarization. 

Information on the ground state geometry of the tran- 
sition metal clusters is usually obtained from the experi- 
ments involving chemical probe methods and photoelec- 
tron spectroscopy, though, such studies for the Co„ clus- 
ters are very limited and not definitive. Reactions of 
Co„ clusters with ammonia and water— indicate icosahe- 
dral structures for the bare and ammoniated clusters in 
the size range n = 50—120 and nonicosahedral packing 
for small (around 19 atoms) Co„ clusters. Although the 
structures of ammoniated Fe„, Co„ and Ni„ clusters in 
the size range of n ~ 19—34 atoms have been found to 
be polyicosahedralfi^ it has been mentioned that the bare 
clusters probably adopt a variety of structures. The pho- 
toionization experiment indicated icosahedral atomic 
shell structures for large Ni„ and Co„ clusters of 50—800 
atoms. However, structures were not well identified for 
small Co„ clusters (n ^50) because atomic sub-shell clos- 
ings in different symmetry based clusters occur in close 
sequences. These experimental results put together indi- 
cate that the icosahedral growth pattern for small sized 
Co„ clusters is less evident. 

Theoretical works on cobalt clusters are limited and 
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the available results are contradictory. Li and Gui^ 
performed first-principles calculation of small Co„ clus- 
ters (4^n^l9) using spin-polarized discrete variational 
method within local density functional theory (DFT). 
However, they had not optimized the structures and con- 
sidered only some special structures with lattice param- 
eters same as the bulk Co. Guevara et alM used an 
unrestricted Hartree-Fock (HF) tight-binding formalism, 
starting from spd-bulk parameterization, but they only 
considered fixed body-centered cubic (bcc) and fee ge- 
ometries for a maximum of 177 atoms without structural 
relaxation. Andriotis and Menorti^ have used a tight- 
binding molecular dynamics scheme to study cobalt clus- 
ters for some selected cluster sizes. Castro et alr^ per- 
formed all-electron density functional calculations using 
both local density and generalized gradient approxima- 
tions. However, the size of the clusters were limited only 
up to 5 atoms. Recently, Lopez et alJ^ studied Co„ 
clusters (4^n^60), where minimization was done using 
an evolutive algorithm based on a many-body Gupta 
potentialii and magnetic properties have been studied 
by a spd tight-binding method. As compared to ab- 
initio methods, the parameterized tight-binding Hamil- 
tonian reduces the computational cost drastically, but 
its main problem is the lack of transferability of its pa- 
rameters. In particular, because of the lack of DFT like 
self-consistency the charge transfer effects are not prop- 
erly accounted for and hence magnetic moment results 
are not fully reliable. 

In this communication, we report the first-principles 
calculation of Co„ clusters (2^n^20). Without any sym- 
metry constraints, we simultaneously relax the geometric 
and magnetic structure to find out the true ground state. 
Our main interest is to study the evolution of structural, 
bonding and magnetic properties as a function of clus- 
ter size. We would also point out how the average bond 
length and coordination number determine the effective 
hybridization and hence the cluster magnetism. 



II. 



COMPUTATIONAL DETAILS 



The calculations are performed using density func- 
tional theory, within the pseudopotential plane wave 
mcthodfi^ We have used projector augmented wave 
(PAW) methodi^i^ and Perdew-Bruke-Ernzerhof (PBE) 
exchange-correlation functional^! for spin-polarized gen- 
eralized gradient approximation (GGA). The 3d and 4s 
electrons are treated as valence electrons and the wave 
functions are expanded in the plane wave basis set with 
the kinetic energy cut-off of 335 eV. Reciprocal space 
integrations are carried out at the F point. Symmetry 
unrestricted geometry and spin optimizations are per- 
formed using conjugate gradient and quasi-Newtonian 
methods until all the force components are less than a 
threshold value of 0.005 eV/A. Simple cubic super-cells 
are used with the periodic boundary conditions, where 
two neighboring clusters are kept separated by at least 



12 A vacuum space. This essentially makes the inter- 
action between the cluster images negligible. For each 
size, several initial geometrical structures have been con- 
sidered. To get the ground state magnetic moment we 
have explicitly considered all possible spin multiplicities 
for each geometrical structure. The binding energy per 
atom is calculated as. 



nE{Co) - E{Con) 



(1) 



where n is the size of the cluster. ^^(Co) and i?(Co„) are 
the total energies of isolated Co-atom and n-atom Co„ 
cluster, respectively. In such a definition, a positive sign 
in El, corresponds to binding. 



III. RESULTS AND DISCUSSIONS 

A. Small Clusters: C02 — Coio 

Both experimental and theoretical predictions of the 
true ground state of the C02 dimcr arc controversial. 
The first experimental estimation of C02 dimer bond 
length and binding energy has been made by mass 
spectroscopy,— which are 2.31 A and 1.72 eV, respec- 
tively. However, more recent collision-induced dissocia- 
tion (CID) experiments^ has estimated an upper bound 
of 1.32 eV to the dimer dissociation energy. The present 
calculation gives dimer binding energy as 1.45 eV/atom 
and a bond length of 1.96 A, which is 78% of the bulk 
hep Co. The Co atoms in dimer have bonding configura- 
tion closer to Sd^As^ than that of the isolated Co atom, 
which is Sd'^As^ and in addition to the highly delocalized 
4s electrons, the more localized id electrons also con- 
tribute strongly to the bonding;^ which consequently, 
produces a shorter bond length for the dimer. Compared 
to the neutral C02 dimer, the experimentall y^^'^^'^^ pre- 
dicted that the bonding in CoJ dimer cation is much 
strong, ~ 2.73 ± 0.27 cV, which is formed by combining 
a neutral Co (3d^4s^) atom with a Co+ (3(i*) cation in 
their respective ground state. Therefore, no promotional 
energy is required to form the cationic dimer. In fact, 
the bonding in CoJ is relatively strong compared with 
other first-row transition metal dimcr cations. On the 
other hand, the formation of neutral C02 dimcr requires 



3d^4s 



3d*4s-'^ promotion^ for both the Co atoms. 



which is 0.42 eV<^ We found that C02 dimer has a total 
magnetic moment of 4 hb, which is also consistent with 
mass spectroscopic measurement^. Our estimates agree 
with the previous first-principles calculations. ^^'^^1'^°'^^ 

In Fig[T] we show the geometrical structures of Co„ 
clusters for the ground state and the first isomer for 
n =2-10 and the calculated binding energy, relative en- 
ergy to the ground state and magnetic moment are given 
in the Table U for the entire size range, n —2 — 20. 

For C03 cluster, we have studied both the linear and 
the triangular structures. An isosceles triangle with to- 
tal magnetic moment 5 /xb is found to be the ground 
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state with binding energy 1.78 eV/atom. Each of the 
two equal sides has length of 2.19 A and other one has 
2.10 A length. Another isosceles triangle with two long 
and one short bond lengths of 2.25 and 2.06 A, respec- 
tively, is found to be nearly degenerate with the ground 
state structure (energy difference is only 3 meV). Ac- 
cording to the present calculation, the linear structure 
lies much higher in energy. The optimal linear structure 
has a total magnetic moment of 7 /^s and lies 0.43 eV 
higher than the ground state. Present result is consistent 
with the spin resonance spectra of C03 in Ar/Kr matrix, 
which indicated a triangular structure with a total mo- 
ment of 5 or 7 to be the ground statci^ The bonding 
in Co;^ cation is stronger than the neutral C03 as it is 
the case for dimeri^ Yoshida et ali^ reported that CoJ 
has a linear structure with bond distance of 2.25—2.50 
A based on their photoelectron spectroscopic study. In 
agreement, previous all-electron (AE) density functional 
calculationi^ predicted an isosceles triangle (2.12, 2.12, 
2.24 A) with a magnetic moment of 1.7 /xs/atom as the 
ground state for C03. This is also consistent with the 
tight-binding study,^^ but they predicted much higher 
bond lengths and magnetic moment. 

Three different initial geometries have been considered 
for C04 cluster: tetrahedral, rectangular and linear. A 
distorted tetrahedron with a total magnetic moment of 
10 fiB appears to be the most stable structure. It has 2.27 
eV/atom energy and has an average bond length of 2.34 
A. Among the six sides of this tetrahedral ground state, 
two pairs have equal lengths of 2.14 A, while the third 
pair is much larger, 2.72 A. There is no experimental re- 
sult available on the structure of neutral tetramer. How- 
ever, Yoshida et alM- predicted a tetrahedral structure 
with a bond length of 2.25 ± 0.2 A as the ground state 
for C04 anion. The initial rectangular structure becomes 
a rhombus after optimization, which also has 10 fiB mag- 
netic moment and lies 0.11 eV higher in energy from the 
ground state. It is the next energetically favorable state, 
which has sides of length 2.14 A and two diagonals of 2.67 
A and 3.35 A. Our results for C04 are consistent with pre- 
vious calculations F^^ii^iii^ii^ii^ Castro et ali^ predicted 
a strong Jahn- Teller distorted tetrahedral ground state 
with bond lengths almost equal to the present values. 
The distorted tetrahedral ground state structure is ac- 
complished by a reduction of some inter-atomic distances 
(and the enlargement of other bonds) until some short 
equilibrium bond lengths result for which there is a more 
effective participation of the (short-range) 3(i-elcctrons 
in the bonds. In fact, in the distorted tetrahedron, there 
are some bonds (these are always on opposite TM— TM 
sides), which have lengths close to that of the dimer. 
These short bonds have high id contributions and they 
are, therefore, the major source of increase of the bond- 
ing in the distorted structure. For both the ground state 
and the first isomer, the magnetic moment is found to be 
2.5 /is/atom. The optimal linear structure is at a much 
higher energy than the ground state. 

We took trigonal bi-pyramid, square pyramid and two 



planar structures: two triangles connected through a ver- 
tex and a pentagon as the initial structures for C05 clus- 
ter. The trigonal bi-pyramid with total magnetic mo- 
ment 13 /is is found to be the most stable structure. This 
structure has 2.55 eV/atom binding energy and 2.34 A 
average bond length. In this ground state (Fig[T]), there 
are two types of bond lengths: all the sides of the upper 
and lower triangular pyramids are of same length and 
are smaller (2.18 A), while those of the interfacing pla- 
nar triangle are much larger, 2.65 A. Another triangular 
bi-pyramid and a square pyramid with an equal mag- 
netic moment of 11 ^b are found to be the degenerate 
first isomer. They lie 125 mcV higher in energy. The 
optimal planar pentagon with 11 ^b magnetic moment 
lies much, 1.04 eV, higher and the double triangle struc- 
ture lies even higher in energy from the ground state. 
Present results are in agreement with the previous AE- 
GGA calculatioufi^ where they predicted the same geo- 
metric structure with 2.28 A average bond length and 2 
/is/atom magnetic moment as the ground state. On the 
other hand, the predictiorti^ of average bond length and 
magnetic moment using Gupta potential is much higher, 
though it predicted the same geometry. 

We have studied the capped trigonal bi-pyramid, oc- 
tahedron and pentagonal pyramid to search the ground 
state for Cog cluster. From now on for the larger clus- 
ters, the planar structures have been discarded by intu- 
ition. An octahedral structure with 14 /i^ total mag- 
netic moment is found to be the ground state. An initial 
capped triangular bi-pyramidal structure relaxes to the 
octahedral ground state. Each side of this octahedral 
ground state is about 2.27 A and has a binding energy 
of 2.93 eV/atom. Another slightly distorted octahedron 
with 12 ^B moment appears as the first isomer. How- 
ever, it is 0.87 eV higher compared to the ground state. 
In the present calculation, the optimal pentagonal pyra- 
mid lies much higher (1.7 eV) in energy compared to 
the ground state and also has 12 magnetic moment. 
Present result is in agreement with previous theoretical 
studiesr^ii^i^i^i and the octahedral structure is gener- 
ally accepted as the most stable structure for Cog clus- 
ter. However, photoelectron spectroscopic study^ pre- 
dicted a pentagonal pyramid with bond distances ~ 2.75 
± 0.1 A to be the most probable structure for Co^ anion 
cluster, i.e., the geometrical structure might strongly be 
correlated with the charged state of the cluster. 

For C07 cluster, we considered capped octahedron, 
pentagonal bi-pyramid and bi-capped triangular bi- 
pyramid. After simultaneous relaxation of both geo- 
metrical and magnetic structure, the capped octahedron 
with 15 ^B magnetic moment appears as the most sta- 
ble structure. This structure has an average bond length 
of 2.29 A and has 2.97 eV/atom binding energy. The 
experimentally measured magnetic moment, 2.36 ± 0.25 
/is/atom^, is little higher than our result. The opti- 
mal pentagonal bi-pyramid, which is a building block of 
icosahedral structure, has a total magnetic moment of 15 
jjLB, which lies 0.19 eV higher in energy from the ground 



state. This is the first isomer and it has an average bond 
length of 2.32 A. The optimal bi-capped triangular bi- 
pyramid has a total magnetic moment of 15 /xb and lies 
0.42 eV higher. However, using Gupta potential, Lopez 
et ali^ predicted a pentagonal bi-pyramidal structure as 
the ground state, and a capped octahcdra as the first 
isomer. 

We have studied three different geometries for Cos 
cluster: bi-cappcd octahedron, capped pentagonal bi- 
pyramid and tri-capped triangular bi-pyramid. The bi- 
capped octahedron with 16 fj,B magnetic moment is found 
to be the most stable structure. This ground state has 
3.07 eV/atom binding energy and an average bond length 
of 2.30 A. The experimentally measured magnetic mo- 
ment, 2.51 ± 0.15 /is/atom;^ is higher than the present 
value. The optimal tri-capped triangular bi-pyramid and 
the optimal capped pentagonal bi-pyramid have an equal 
magnetic moment of 16 but lie 0.4 and 0.48 eV higher 
in energy respectively. They are the first and second iso- 
mers (see FigH]). 

For the Cog cluster, we considered tri-capped octahe- 
dron and bi-capped pentagonal bi-pyramid as initial con- 
figurations. A distorted tri-capped octahedron is found 
to be the most stable structure with 3.14 eV/atom bind- 
ing energy. This ground state structure has 17 jU_b to- 
tal magnetic moment, which is smaller than the ex- 
perimental value of 2.38 ± 0.11 /jb /sAom^ The opti- 
mal bi-capped pentagonal bi-pyramid has 17 jiB mag- 
netic moment and lies 0.53 eV higher in energy. This 
structure is found to the first isomer. The present re- 
sults are not in good agreement with the semi-empirical 
predictions where they predicted the bi-capped pen- 
tagonal bi-pyramid as the ground state with a relatively 
high magnetic moment. 

Different tri-capped pentagonal bi-pyramid (TCPBP) 
structures along with different tetra-capped octahedral 
structures were taken as initial structures for Coio clus- 
ter. A TCPBP structure with 18 hb total magnetic mo- 
ment is found to be the ground state. The calculated 
magnetic moment in the ground state is smaller as com- 
pared to the neighboring sizes, which is indeed the case 
in experiment (cf. FigE] and will be discussed later). 
This is because of the fact that TCPBP is an icosahe- 
dral fragment based on pentagonal bi-pyramid. This is 
different from the structural growth seen for Cog— Cog 
clusters, where the ground state structures are all octa- 
hedral based. For this TCPBP ground state, average co- 
ordination and average bond lengths are slightly higher 
and the competing effect of these two makes the mag- 
netic moment smaller than its neighboring clusters. An- 
other TCPBP with total magnetic moment 20 fis lies 
0.08 eV higher in energy compared to ground state is 
found to be the first isomer. The experimental magnetic 
moment (2.07 ± 0.10 /is/atom^i) is larger compared to 
the predicted ground state but very close to the first iso- 
mer, which is energetically very close to the ground state. 
Lopez et alt^ predicted the same TCPBP structure, but 
with much higher, 2.45 /XB/atom, magnetic moment as 
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FIG. 1: (Color online) Equilibrium geometries of the two ener- 
getically lowest isomers of Con clusters for n = 2—10 (except 
n = 8, for which three consecutive isomers have been shown). 
Numbers in the parenthesis represent number of atoms in the 
cluster, relative energy to the respective ground state and to- 
tal magnetic moment, respectively. 



the ground state, whereas Guevara et alw^ predicted a 
fee ground state with comparable, 2 [Ib /atom, magnetic 
moment. 
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TABLE I: Binding energy, relative energy to the theoretically computed minimum energy state {AE — E — En-iin) and magnetic 
moment for Co„ [n =2—20) clusters. Recent SG experimental results^ i^i'^^i'^^ of the magnetic moment are shown for comparison. 



Cluster 




AE 


Magnetic Moment 


Cluster 


Et 


AE 




Magnetic Moment 




(eV/atom) 


(eV) 


Theory 


(/is /atom) 
SG Exp. 




(eV/atom) 


(eV) 


Theory 


(/is /atom) 
SG Exp. 


C02 


1.452 


0.000 


2.00 




Cois 


3.397 


0.000 


2.07 


2.38±0.03" (2.09±0.04)'' 


C03 


1.783 


0.000 


1.67 


_ 




3.393 


0.046 


1.93 






1.783 


0.003 


2.33 






3.388 


0.125 


2.20 




C04 


2.274 


0.000 


2.50 







3.385 


0.169 


1.80 






2.248 


0.106 


2.50 






3.385 


0.171 


1.93 




Cos 


2.553 


0.000 


2.60 


_ 


C016 


3.458 


0.000 


2.13 


2.53±0.04" (2.32±0.01)'' 




2.528 


0.125 


2.20 






3.458 


0.005 


2.00 




C06 


2.929 


0.000 


2.33 







3.445 


0.208 


1.88 






2.784 


0.869 


2.00 






3.439 


0.308 


2.25 




C07 


2.971 


0.000 


2.14 


2.36±0.25'' 




3.438 


0.319 


1.88 






2.944 


0.192 


2.14 




C017 


3.514 


0.000 


2.06 


2.24±0.04" (2.19±0.02)'' 


Cos 


3.074 


0.000 


2.00 


2.51±0.15" 




3.506 


0.123 


2.18 






3.024 


0.400 


2.00 






3.504 


0.167 


1.94 






3.013 


0.484 


2.00 






3.490 


0.407 


1.82 




Cog 


3.143 


0.000 


1.89 


2.38±0.11'' 




3.466 


0.812 


2.06 






3.084 


0.527 


1.89 




Coig 


3.555 


0.000 


2.00 


2.07±0.04" (2.37±0.07)'' 


Colo 


3.137 


0.000 


1.80 


2.07±0.10" 




3.554 


0.024 


2.11 






3.128 


0.085 


2.00 






3.544 


0.194 


1.89 




Con 


3.205 


0.000 


1.91 


2.42±0.09" 




3.523 


0.571 


2.00 






3.203 


0.016 


1.91 




C019 


3.607 


0.000 


2.05 


i.ziztO.Oo (2.48±0.04j 


C012 


3.252 


0.000 


2.00 


2.26±0.08" (2.21±0.01)' 




3.597 


0.174 


1.95 






3.243 


0.103 


1.89 






3.581 


0.478 


1.84 




C0l3 


3.279 


0.000 


1.92 


2.30±0.07" (2.00±0.06)' 




3.559 


0.901 


1.74 






3.268 


0.140 


2.08 






3.546 


1.158 


2.16 






3.266 


0.167 


2.38 






3.542 


1.220 


1.95 




Coi4 


3.323 


0.000 


2.00 


2.29±0.06" (2.11±0.02)' 


C020 


3.620 


0.000 


2.00 


2.04±0.05" (2.36±0.02)'' 




3.322 


0.004 


2.00 






3.607 


0.262 


1.90 






3.322 


0.005 


1.71 






3.588 


0.634 


1.80 






3.322 


0.007 


1.86 






3.576 


0.891 


2.10 






3.320 


0.008 


2.14 






3.565 


1.103 


1.90 





"From Knickclboin (Rcf.8 and Rcf.35) 
''From Xu et al. (Ref.7 and Ref.36) 



B. Intermediate size clusters: Con C020 

With the increase in the number of atoms in the clus- 
ter, the determination of the ground state becomes a dif- 
ficult task as the number of local minima in the potential 
energy surface increases very rapidly with the number of 
atoms in the cluster. In the search for the ground state 
structures, we have considered several possible geomet- 
rical structures as initial guess and wc relax all of them 
for all possible spin multiplicities for each ri-atom cluster. 
The predicted minimum energy structure along with an 
higher energetic isomer for each cluster in the size range 
n =11—20 have been shown in Figl2l 

The 13- atom hexagonal close packed structure consists 
of a hexagonal ring around a central atom and two tri- 
angular planes above and below it: 3,7,3 stacking, and 
the 13-atom close packed icosahcdral structure has two 
pentagonal rings, two apex atoms and a central atom: 
1,5,1,5,1 stacking. The initial structures for both the 
Coii and C012 clusters, have been derived from these 13- 



atom close packed structures by removing 1 or 2 atoms, 
respectively. As initial configurations for Con, we con- 
sidered three hep derived structures (4,7; 3,5,3 and 1,7,3 
stacking) and two icosahedral derived structures (5,1,5 
and 1,5,1,4 stacking). On relaxation, all the structures 
distorted heavily. The distortion of the 4,7 hexagonal 
structure is such that one atom from the hexagonal ring 
comes out of the plane (Fig[2]) and is the most stable 
configuration with binding energy 3.20 eV/atom. This 
structure has a total magnetic moment of 21 fj.B, which is 
considerably smaller than that of the experimental value, 
2.42 ± 0.09 /is/atom^ The initial 1,5,1,4 icosahedral 
structure is found to be the next isomer after relaxation. 
This structure lies just 16 meV higher in energy from the 
ground state and also has 21 magnetic moment. 

For C012 cluster, we have considered two hep struc- 
tures: one with 3,6,3 stacking (without the central atom) 
and the other with 2,7,3 stacking (with the central atom) 
and two icosahedral structures: a closed ieosahedra with- 
out one apex atom and a closed ieosahedra without the 
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FIG. 2: (Color online) The ground state and a higher energy structure of Co,i clusters (n =11 — 20). For the size range 
n =13—20, we show the optimal icosahedra-derived structures for comparison. The first entry within the parenthesis, n.k, 
indicates that the structure corresponds to the fc-th isomer of Co„ cluster, second and third entries give the relative energy to 
the ground state and magnetic moment per atom, respectively for the fc-th isomer, = corresponds to the ground sate. 



central atom. After relaxation the initial hep structures 
undergo a considerable rearrangement, which consists of 
a plane of 7 atoms coupled with another 5-atoms plane 
(see Fig 12). This structure having a total magnetic mo- 
ment of 24 /is is the lowest energy state, which has 
3.25 eV/atom binding energy. A structural rearrange- 
ment has also been seen for 1,5,5,1 icosahcdral structure, 
which looks like a hep fragment of 3,5,4 stacking after 



relaxation. This structure has a total magnetic moment 
of 22 i^lb and has 3.24 eV/atom binding energy. The 
icosahedral structure with 1,5,1,5 stacking does not lead 
to such rearrangement after relaxation. However, the 
relaxed structure lies much higher in energy compared 
to the minimum energy state. The calculated magnetic 
moment in the lowest energy structure, agrees with the 
recent SG expcrimentsii^ and the previous theoretical 



7 



calculations 

We considered the icosahedral, hep, cuboctahedral and 
fee structures as initial structures for C013 cluster. A 
distorted hexagonal structure is found to be the mini- 
mum energy state and it looks like the most stable Con 
structure along with two additional capped atoms (see 
Figl2]). This structure has a total magnetic moment 
of 25 fiB and has 3.28 eV/atom energy. Knickelbein 
found the experimental magnetic moment to be 2.30 ± 
0.07 ^s/atom,— which is slightly higher than the present 
value, 1.92 yUs/atom. However, this calculated value is in 
good agreement with another recent SG experiment by 
Xu et al.-^^ which predicted 2.00±0.06 /j,s/atom mo- 
ment. Another distorted hep structure with total mag- 
netic moment 27 /is is found to be the first isomer, which 
lies 0.14 eV higher in energy from the ground state. The 
optimal icosahedral structure has a comparatively large 
magnetic moment, 31 /zb, and it lies 0.17 eV higher in 
energy being the second isomer. The optimal fee and 
cuboctahedral structures are much higher in energy with 
respect to the minimum energy state. The present pre- 
diction of hep structure as the minimum energy state is in 
agreement with the previous tight-binding predictions^ 
However, there arc some reports,—!^ which favor the 
icosahedral structure as the ground state. 

For C014 cluster, the trial structures are complete 
icosahedra with a single atom capping and a hexagonal 
structure with 3,7,4 stacking. The optimal icosahedral 
and the optimal hexagonal structures are found to be 
degenerate. The energy separation is only 4 meV. Both 
of these structures have equal (28 ijlb) magnetic moment 
(see Fig 12]). We also found several isomers which lie very 
close to these structures: An icosahedral structure (24 
/Ub), a hexagonal structure (26 /is) and another icosahe- 
dra (30 ^b) lie only 5, 7 and 8 meV above the ground 
state, respectively (see Table The very recent SG 
experimental predictions range magnetic moment from 
2.11± 0.02 ^B /atomize to 2.29 ± 0.06 ^slatouS^ for 
C014 cluster, which is in good agreement with the present 
result. 

The C015 trial structure with lowest energy is the 
hexagonal structure with 4,7,4 atomic staking. This 
structure has 31 ^b magnetic moment, which is in agree- 
ment with the SG experiment of Xu et ai, which is 
2.09±0.04 /is/atom.— However, Knickelbein predicted 
a larger value.— The other hexagonal structures with 
total magnetic moments 29, 33 and 27 hb lie ~ 0.05, 
0.12 and 0.17 eV higher than the minimum energy 
state, respectively. The optimal icosahedral structure of 
1,5,1,5,1,2 stacking (with 29 /is magnetic moment) is the 
fourth isomer, which lies 0.17 eV above the lowest energy 
state. 

The same kind of structural growth is observed in the 
case of C016 cluster. The hexagonal structure with a 
total magnetic moment of 34 fiB is found to be the low- 
est in energy. This structure has 4,7,5 stacking and 3.46 
eV/atom binding energy. This structure is nearly degen- 
erate (5 meV lower) with another hexagonal structure, 



which has 32 ^b magnetic moment. The next two iso- 
mers are also of same hexagonal motif, which have 30 
and 36 /is magnetic moment and they lie 0.21 and 0.31 
eV higher in energy with respect to the lowest energy 
state, respectively. The optimal icosahedral structures 
with 1,5,1,5,1,3 stacking and with magnetic moments 30 
and 32 fiB lie 0.32 and 0.36 eV higher, respectively. The 
optimal hexagonal and icosahedral structures are shown 
in FiglH The other icosahedral structure with 5,1,5,4 
stacking is found to be much higher in energy. 

The hexagonal structure with total magnetic moment 
35 /is is the lowest energy state for C017 cluster. This 
structure has 5,7,5 stacking and 3.5 cV/atom binding en- 
ergy. The calculated magnetic moment, 2.06 /^s/atom, 
is slightly smaller than that of predicted by both the re- 
cent experimentsii^ (sec Table|I|. The next three isomers 
also have hexagonal symmetry. They have 37, 33 and 31 
fiB magnetic moment and they lie 0.12, 0.17 and 0.41 eV 
higher than the lowest energy state, respectively. The op- 
timal icosahedral structure (Fig[2]) has 5,1,5,1,5 stacking 
and lies much higher in energy. 

A 6,7,5-hcp trail structure (FiglJ), which has a total 
magnetic moment of 36 fiB is found to be the lowest 
energy state for Cois cluster. This structure has a bind- 
ing energy of 3.61 eV/atom. The magnetic moment is 
in agreement with Knickelbein, 2.07±0.04 /is /atom<^i2^ 
However, Xu et al. predicted a higher value, 2.37±0.07 
Hb/ atomiii^^ Another two hep structures with total mag- 
netic moments 38 and 34 hb lie 0.02 and 0.19 eV higher in 
energy, respectively and are the first and second isomers. 
The optimal icosahedral structure has a total magnetic 
moment of 36 jiB, which lies much higher (0.57 cV) in 
energy and is the third isomer. 

For C019 cluster, we investigated a double icosahe- 
dral structure, a hep structure with 6,7,6 stacking and 
a cuboctahedral structure. The hep structure (FigE]) 
with 39 /is magnetic moment appears as the most stable 
structure among all the trial structures. The calculated 
magnetic moment, 2.05 /is/atom, is closer to the value 
of Knickelbein, 2.21±0.03 ^s/atom^i^ than that of the 
value predicted by Xu et al, 2.48±0.04 /is/atom,!^ The 
next four isomers are also found to be of same hep pack- 
ing. These isomers with total magnetic moments 37, 35, 
33 and 41 lie 0.17, 0.48, 0.90 and 1.16 eV higher than 
the minimum energy state, respectively. On the other 
hand, the optimal icosahedral structure (FiglJ) has 37 
fiB magnetic moment and lies 1.22 eV higher from the 
ground state. The fee and hep fragments have also been 
proposed as ground state structures in the previous theo- 
retical calculationsJ^ii^ Our results arc in according with 
the predicted hep structure. Also some calculation o^^i^^ 
predicted icosahedral ground state for Coig. 

Among all the considered structures the capped 19- 
atom hexagonal structure (Fig[21) is found to be the low- 
est energy state for C020 cluster. The calculated mag- 
netic moment is found to be 2 /^s / atom for this struc- 
ture, which is in agreement with the value predicted by 
Knickelbein^ (2.04±0.05 //s/atom). However, the mo- 
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ment predicted by Xu et al. is much higher (2.36±0.02 
/is/aton>i). Similar to what we have seen for C015— C019 
clusters, the next few isomers are also of hexagonal mo- 
tif. The hep structures with total magnetic moments 38, 
36, and 42 ^b, which are 0.26, 0.63 and 0.89 eV higher 
are found to be the first, second and third isomers, re- 
spectively. The optimal icosahedral structure has a total 
magnetic moment of 38 and appears as the forth iso- 
mer (Fig[2]). However, this structure lies much higher in 
energy. 



C. Binding energy, stability and dissociation 

energy 

Calculated binding energies are plotted in FigOfor the 
ground states of Co„ clusters in the size range n =2—20. 
Since the coordination number increases with the num- 
ber of atoms in the cluster, the binding energy increases 
monotonically. The binding energy of the largest cluster 
studied here (C020) is 3.62 eV/atom, which is about 82% 
of the experimental bulk value, 4.4 eV/atom^^l for hep 
Co. Upon extrapolation of the linear fit of the binding 
energy per atom data to n~^/'^ — > (Fig[3Ja)), we can 
estimate the binding energy of the infinitely large clus- 
ter. This is found to be 5.0 cV/atom, which is larger 
than the experimental value for hep bulk Co. However, 
within the same level of theory we found the hep bulk co- 
hesive energy to be 5.11 eV/atom, which is close to the 
extrapolated value but again larger than the experimen- 
tal value. This overestimation is consistent with the DFT 
calculationii Calculated binding energies for the optimal 
hexagonal and optimal icosahedral structures (see Fig[2] 
for the optimal geometries) and the corresponding energy 
difference (i?diff ) between them are plotted in FigH] for 
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FIG. 3: Plot of binding energy per atom as a function of 
cluster size n for the ground state structures of the entire size 
range 2^ n ^20. (a) Plot of binding energy per atom (B.E) 
as a function of n~^^^ for the clusters Co„, 6^ n ^20 and a 
linear fit (-3.90 n~^''^ + 5.00) to the data, (b) Plot of second 
difference in total energy (A2-B), which represents the relative 
stability. 
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FIG. 4: (color online) Plot of binding energy per atom (_£(,) for 
optimal hexagonal and optimal icosahedral structures (left) 
and plot of total energy difference between these optimal 
structures, iJdiff = —[£■( hexagonal) — iJ (icosahedral)], (right) 
for the size range n =13—20. -Ediff increases with cluster size. 



the size range n =13—20. The optimal hexagonal struc- 
tures are always found to be the ground state for this 
size range except for Co 14, where the optimal hexagonal 
and icosahedral structures are found to be degenerate. 
Moreover, in this size range, next few isomers are also of 
hexagonal motif and the optimal icosahedral structures 
appear as higher energy (third, fourth or fifth) isomers for 
n =15—20 (see FiglS]). The energy difference between the 
hexagonal ground state and optimal icosahedral struc- 
tures increases with increasing cluster size making icosa- 
hedral structures more and more unfavorable. We plot 
the energy variation as a function of cluster magnetic 
moment for icosahedral and hexagonal C013, C015, C017, 
and Coig clusters in FigEl Both the structures show sim- 
ilar qualitative behavior for all the clusters and they have 
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FIG. 5: (color online) Plot of total energy as a function of 
magnetic moment for C013 (a), C015 (b), C017 (c) and C019 
(d) clusters. The • (■) represents icosahedral (hexagonal) 
structure. Insets represent magnification around the minima. 
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hexagonal minima around 2 /is/atom moment. 

We calculate the second difference in the total energy: 

A2E{n) ^ E{n+l) + E{n-l)~2E{n), (2) 

where E{n) represents the total energy of an n— atom 
cluster. Calculated A2E has been plotted in Fig[3Jb), 
where we see the peaks at n = 6 and 9, i.e., the clusters 
with 6 and 9 atoms are particularly more stable than 
their neighboring clusters. The stable structure for Cog 
is a octahedron and for C09, it is a distorted tri-capped 
octahedron. The CID experiment^ has also been indi- 
cated a maximum at n = 6 in the measured dissociation 
energy, which indicates a higher stability of the hexamer. 
The extra stability of hexamer indicates that the octahe- 
dral structure can act as a building block for larger size 
clusters and, indeed, for C015— C020 clusters, we have 
found a distinct hexagonal growth pattern and an octa- 
hedron is just a fragment of a hexagonal structure. The 
calculated stability (Figl3] (b)) shows minima at n = 3, 
5, 7, 10 and 14, which are related to their weak bonding. 

This can be further demonstrated by studying the dis- 
sociation energies as an 71-atom cluster fragments into m 
and (n — m)-atom clusters. The m-channel dissociation 
energy can be calculated as. 



D,n{n) = E{m) + E{n - m) - E{n), 



(3) 



where E(n), E{m) and E(n — m) are the total ener- 
gies of n, m and {11 — m) atom clusters, respectively. 
We have plotted the calculated single channel {Di) and 
dimer channel (-D2) dissociation energies in Fig|S]and Di 
is compared with the CID experiment by Hales et aLM- 
However, they have estimated this dissociation energy 
through an indirect method: actually, they measured the 
single channel dissociation energy of Co+ cation cluster 
and derived the same for the neutral one by using the 
ionization energies (IE) of the neutral clusters measured 
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FIG. 6: (Color online) Plot of single channel, Di and dimer 
channel, D2, dissociation energies as a function of cluster size 
n for the GS configurations. We compare our calculated single 
channel dissociation with the CID experimental result.— 



TABLE II: Theoretically calculated single channel bond dis- 
sociation energies(BDE) compared with experimentally mea- 
sured values^^ for Co„ (n =2—20). Experimental uncertain- 
ties are within parentheses. 



Co„ 


BDE(eV) 


Co„ 


BDE(cV) 




Theory CID Exp. 




Theory CID Exp. 


2 


2.90 


<1.32 


12 


3.77 


3.41(0.28) 


3 


2.45 


>1.45 


13 


3.61 


3.68(0.31) 


4 


3.75 


2.41(0.21) 


14 


3.89 


3.12(0.36) 


5 


3.67 


2.84(0.28) 


15 


4.43 


3.84(0.39) 


6 


4.81 


3.31(0.29) 


16 


4.38 


3.82(0.39) 


7 


3.22 


2.65(0.24) 


17 


4.40 


3.44(0.44) 


8 


3.79 


2.93(0.22) 


18 


4.26 


3.84(0.56) 


9 


3.70 


2.89(0.17) 


19 


4.54 




10 


3.08 


3.05(0.16) 


20 


3.88 




11 


3.89 


3.12(0.26) 









by Yang and Knickelbein,— and Parks et al.^ i.e., 
i^fP = i?i(Co+) + /^(Co„) - /S(Co„_i). 



(4) 



The calculated single channel dissociation energy, Z?i, 
shows a high peak at n =6 and dips at n = 5, 7 and 10, 
which are consistent with our stability analysis. However, 
we do not find any dip in the calculated dissociation en- 
ergy at n =14, as has been seen in the CID experiment. 
Generally, the single channel dissociation energy is the 
most favorable except for n =4, where the dimer disso- 
ciation (C04 C02 -f C02) is more favorable than the 
single channel (C04 C03 4- Co) dissociation. Table 
II shows the theoretically computed single channel bond 
dissociation energy compared to the experimentally mea- 
sured values^ for the entire range of clusters having sizes 
2 to 20. 

To understand the optimized structures, we calculated 
the average bond lengths and average coordination num- 
ber for the ground state geometries and plotted them 
in FiglTlJa) and FiglJUb), respectively, as a function of 
cluster size. These two quantities are closely related to 
the structure of the cluster. We define the average bond 
length as (r) = ^ Si> j ''u ' where is the bond dis- 
tance between the j-th and i-th atoms, and Ub is the 
number of such bonds. Here we consider that two atoms 
are bonded if their inter atomic distance is within 2.91 
A, which is around the average of the first (2.51 A) and 
second (3.54 A) nearest-neighbor distances in bulk Co. 
The average coordination number in a cluster is defined 
as (ric) = i where is the number of neighbors 

within the chosen cut-off of the fc-th atom in the cluster 
of n atoms. The convergence of the average bond length 
to the bulk value (2.51 A) is much faster than the con- 
vergence of average coordination, which is far below the 
bulk value (12 for hep Co). Dips at n = 6 and 9 in Fig. 
[Tlja) indicate that the atoms in these clusters are closely 
spaced and strongly bonded compared to the neighbors, 
and therefore are more stable than the neighboring struc- 
tures. While the peaks at n = 5, 10 and 14 in Fig[7|^a) 
and at n =10 and 14 in FiglTl^b) indicate that atoms in 
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FIG. 7: Plot of (a) average bond length (r) and (b) average 
coordination (nc), as function of cluster size ?i for the ground 
state geometries. Dots and squares represent the data points, 
while solid line through them is guide to eye. 

these clusters are far apart and slightly more coordinated 
then their neighbors)^ which results in a weak bonding 
in these clusters compared to their neighbors. 



D. Magnetic moment 

The calculated magnetic moments are plotted in Fig[5] 
as a function of cluster size (n). The Co— Co interaction 
is always ferromagnetic for the entire size range studied, 
as it is for hep bulk Co. However, the magnetic moment 
(2—2.5 /xs/atom) is larger than the hep bulk value, 1.72 
Hb/ atomj^ This enhancement in moment for a few atom 
cluster can readily be understood from the more local- 
ized d-electrons resulting from the decrease in effective 
hybridization. The calculated magnetic moments are in 
fair agreement with the very recent SG experiments by 
Xu et al^ and Knickelbein.— Fig[8] shows a qualitative 
agreement between the calculated and the experimental 
values^i^i^i^^ though the calculated moments are always 
underestimated systematically. However, calculated mo- 
ments are close to that of predicted by Xu et al^^ for 
the size range n =13—17 and in the size range n =18—20 
they are close to the values predicted by Knickelbein 
The underestimation of calculated moment may be due to 
the fact that we did not include spin-orbit interaction in 
the present calculation. Moreover, one should remember 
that the magnetic moments in a magnetic deflection mea- 
surement are always derived assuming a model, which 
may influence the outcome. For example, Knickelbein^ 
used either superparamagnetic or locked moment model, 
whether Xu et al^ assumed an adiabatic magnetization 
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FIG. 8: (Color online) Size dependent variation of magnetic 
moment of the corresponding ground states. Calculated mag- 
netic moments of the first isomers have also been shown. Cal- 
culated magnetic moments are compared with experimental 
results (Ref.8 and Ref.7). 



model to derive the moments experimentally for cobalt 
clusters. It is important to note that both experiments 
show same size evolution in general but there are some 
systematic differences. However, this is not due to the 
adoption of different models to calculate the magnetic 
moment as both the models resemble with the same Curie 
law for magnetization,i^ but may be due to differing iso- 
mer distribution in the SG beam. 

The magnetic moment is strongly correlated with the 
effective hybridization, which is closely related to the av- 
erage bond length (r) and the average coordination num- 
ber (ric). As {uc) decreases the magnetic moment should 
increase through the decrease in effective hybridization. 
On the other hand, the dependency of magnetic moment 
on (r) is directly proportional: a decrease in (r) results in 
decrease in magnetic moment through the enhancement 
in effective bonding. FiglTja) and FiglTl^b) show that as 
we go from n — A to n = 10, both (r) and (ric) increase, 
whereas Fig|S]shows that the magnetic moment per atom 
decreases. Therefore, between these two competing con- 
tributions ((r) and {ric)) to the magnetic moment, the 
average coordination number dominates over the average 
bond length in the size range n =4—10. 

In the intermediate size range, n = 11—20, the vari- 
ation of (r) (Figl?^) and (n^) (FiglTjD) is much slower 
with n, and therefore, the magnetic moment per atom 
does not vary rapidly. It is around 2 fis/ atom for all 
the clusters in this size range. So, in this size range, 
it is hard to predict the dominant parameter for mag- 
netism. To illustrate the effect of (r) and (uc) on the 
magnetism in the intermediate size range we compare 
these two quantities for the optimal hep and icosahedral 
structures (sec insets (a) and (b) of FiglS]). It is seen that 
for a hep structure, both the (r) and (ric) arc smaller than 
those of corresponding icosahedral structure for a partic- 
ular n-atom cluster. In addition the magnetic moments 
of optimal hep clusters are always larger than or equal 
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FIG. 9: (Color online) Comparison of magnetic moment be- 
tween optimized hep and optimized icosahedral structures for 
Co„ in the size range n =15—20. The filled squares and filled 
circles correspond to results for hep and icosahedral struc- 
tures, respectively. The insets show the corresponding com- 
parisons for (a) average bond length and (b) average coordi- 
nation number. 



n =2—20. In the intermediate size range, the clusters 
adopt hep structural packing, which is different from 
the trend observed for the other 3d transition metal 
clusters .i^i*^ In the size range n =15—20, the energy dif- 
ference between the hexagonal minimum energy states 
and optimal icosahedral structures increases with cluster 
size making the icosahedral structures more and more un- 
favorable with increasing size. The calculated magnetic 
moments are in good agreement with both the recent 
SG experiments. It is found that the effect of average 
coordination number always dominates over the average 
bond length to determine the effective hybridization and 
therefore, the magnetic moment of the clusters. 



to that of the corresponding optimal icosahedral clusters 
(see Figl9]), which again demonstrates that in this inter- 
mediate size range also the coordination dominates over 
the average bond length in deciding magnetism. 



IV. SUMMARY AND CONCLUSIONS 

We have systematically studied the structure, bond- 
ing and magnetism in Co„ clusters in the size range 
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